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Fluctuation measurements of event-wise mean transverse momentum (pt) for p-p and Pb-Pb colli¬ 
sions at the large hadron collider (LHC) have been reported recently. In that study it was concluded 
that the strength of “nonstatistical” ( p t ) fluctuations decreases with increasing particle multiplicity 
n c h (or A-A centrality) and is nearly independent of collision energy over a large interval. Among 
several potential mechanisms for those trends onset of thermalization and collectivity are men¬ 
tioned. The LHC analysis employed one fluctuation measure selected from several possibilities. 

An alternative fluctuation measure reveals strong increase of pt fluctuations with n c h (or A-A cen¬ 
trality) and collision energy, consistent with previous measurements at the relativistic heavy ion 
collider (RHIC). The pt fluctuation data for LHC p-p collisions can be described accurately by a 
two-component (soft+hard) model (TCM) in which the hard component represents dijet produc¬ 
tion. The data for Pb-Pb collisions are described accurately by a TCM reference for more-peripheral 
collisions (suggesting transparent collisions), but the data deviate quantitatively from the reference 
for more-central collisions suggesting modification of jet formation. Overall fluctuation data trends 
suggest that minimum-bias jets (minijets) dominate pt fluctuations at both the LHC and RHIC. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 


I. INTRODUCTION 

Fluctuation measurements of event-wise mean trans¬ 
verse momentum, denoted in this study by ( p t ), vs charge 
multiplicity n c h in p-p and Pb-Pb collisions at the large 
hadron collider (LHC) have been reported recently PQ. A 
principal motivation for such measurements is the search 
for evidence of the phase transition between a quark- 
gluon plasma (QGP) and a hadronic medium (hadron 
gas) in the form of excess fluctuations of a thermody¬ 
namic quantity, specifically event-wise (pt) as a proxy 
for a local temperature [2j. The (pt) fluctuation mea¬ 
sure chosen from among several candidates is apparently 
based on the assumed thermodynamic context. The re¬ 
ported systematic behavior includes negligible energy de¬ 
pendence over a large energy interval and general de¬ 
crease of fluctuation “strength” with increasing event 
multiplicity or centrality. Those results conflict sharply 
with previous relativistic heavy ion collider (RHIC) re¬ 
sults employing alternative fluctuation measures in which 
p t fluctuation amplitudes increase strongly with Au-Au 
collision centrality and with collision energy EH5]- 

It is reasonable to assume direct connections between 
underlying collision mechanisms and final-state collision 
structure in the form of yields, spectra and correlations. 
But how the structure should be characterized statis¬ 
tically and how results should be interpreted in terms 
of physical mechanisms is in question. The overarching 
goal should be consistent and interpretable descriptions 
of particle densities on a multidimensional momentum 
space varying event-wise over collision events and extrac¬ 
tion of all information residing in such data distributions. 

The present study focuses on fluctuation measurement. 
Fluctuations and angular correlations are intimately re¬ 
lated m©. The scale (bin size) dependence of fluctu¬ 


ations corresponds to the space distribution of angular 
correlations [7]. The more structured a multiparticle an¬ 
gular distribution the larger the event-wise fluctuations. 
Fluctuation measures should then be compatible with 
correlation measures, providing useful design constraints. 

Just as for Ref. [IJ some previous fluctuation measure¬ 
ments were motivated by a search for excess or critical 
fluctuations near a QCD phase boundary [SHID]- Con¬ 
jectured possibilities included (a) QGP formed in some 
special events comprising a small fraction of an event en¬ 
semble or (b) QGP formed in most events of some event 
class resulting in a general fluctuation excess relative to 
some statistical reference. While (a) was not observed 
and evidence for (b) was marginal at the SPS [9], excess 
mean-p t fluctuations in at least most events were clearly 
apparent in the first Au-Au data from the RHIC m- 
Subsequent measurements of pt fluctuations, particu¬ 
larly 2D bin-size dependence, led to a surprising result: 
p t angular correlations inferred from inversion of p t fluc¬ 
tuation bin-size dependence revealed that dijet produc¬ 
tion is the dominant source of pt fluctuations at RHIC 
energies |4j. The collision-energy dependence of p t fluc¬ 
tuations and inferred angular correlations is also consis¬ 
tent with QCD dijet production via parton scattering [3]. 
Those results seemed consistent with number (distin¬ 
guished from pt) angular correlations that similarly in¬ 
dicated a dominant role for dijet production dunl¬ 
in contrast, fluctuation measurements employing al¬ 
ternative statistical measures (per-pair rather than per- 
particle, those terms defined below) applied to RHIC 
data suggested a very different scenario: p t fluctua¬ 
tions/correlations, in ratio with ensemble-mean pt , de¬ 
crease strongly with increasing A-A centrality (suggest¬ 
ing increased thermalization) and are essentially indepen¬ 
dent of collision energy m The recent analysis of LHC 
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p-p and Pb-Pb data [T] seems to confirm those results. 

The present study is the followup to a recent analy¬ 
sis of ensemble-mean p t (denoted here by p t ) systemat- 
ics m , derived from LHC p-p, p -Pb and Pb-Pb data m 
wherein the p t data for three collision systems and several 
energies are accurately described by a two-component 
(soft+hard) model (TCM) featuring dijet production as 
the hard component predicted quantitatively by pertur¬ 
bative QCD (pQCD). In contrast, theory Monte Carlos 
based on final-state hadron or parton rescattering and/or 
transverse flows US] fail to describe the p t (n c h ) data. 

Resolution of the pt fluctuation dichotomy requires de¬ 
tailed analysis of methods and interpretations. I first re¬ 
view the definitions and properties of several statistical 
measures applied to event-wise p t fluctuations, with ref¬ 
erence to measure design criteria. I present results from 
the LHC fluctuation analysis reported in Ref. [T] in alter¬ 
native plotting formats, discussing the apparent physical 
implications of the various data trends. Finally, I relate 
the LHC results to RHIC results derived from p t fluctua¬ 
tion inversion to pt angular correlations and to measured 
dijet systematics. The evidence suggests that an existing 
comprehensive TCM scenario describing dijet manifesta¬ 
tions in yields, spectra and correlations also describes pt 
fluctuation data accurately and supports the conclusion 
that dijet production plays a central role in all aspects of 
hadron production in high-energy nuclear collisions. 

This article is arranged as follows: Section [TT] sum¬ 
marizes a recent study of (p t ) fluctuations at the LHC. 
Section m introduces various methods for p t fluctuation 
measurement. Section m defines several basic models 
for p-p and A-A collisions. Section |V| defines a p-p TCM 
for ensemble-mean p t . Section m summarizes LHC p-p 
p t fluctuation data. Section |VII presents LHC Pb-Pb p t 
fluctuation data compared to RHIC Au-Au data. Sec¬ 
tion m relates pt fluctuation measurements to pt an¬ 
gular correlation measurements including dijet-related 
structure. Sections |IX| and |X| present discussion and sum¬ 
mary, and App. [A] presents a pt TCM for A-A collisions. 


II. LHC (p t ) FLUCTUATION ANALYSIS 

Reference |T] (ALICE collaboration) reports measure¬ 
ments of ( p t ) fluctuations vs n c h in p-p and Pb-Pb colli¬ 
sions at the LHC. Here I summarize conclusions of that 
study and in following sections consider details of p t fluc¬ 
tuation trends. In Ref. [T] event-wise mean-p t is denoted 
by MEbE(PT) whereas in this study the symbol is ( p t )■ 
In the context of a temperature narrative the primary 
motivation for ( p t ) fluctuation analysis is a search for crit¬ 
ical fluctuations of (local) temperature associated with a 
QCD phase boundary. A QCD phase transition or crit¬ 
ical point may “go along with” critical fluctuations in 
a thermodynamic quantity such as temperature repre¬ 
sented (within an assumed theoretical context) by ( pt ). 
The adopted dimensionless ratio \/Cm/pt,m is said to 
“quantify the strength of the non-statistical fluctuations 


in units of the [ensemble] average transverse momentum 
M(pT)m [= Pt,m] in the multiplicity class m.” 

The basic fluctuation measure is the ratio C = 
B/ n c h{n c h — 1), where I? is a variance difference defined 
below. The p t acceptance is p t G [0.15,2] GeV/c, where 
the lower limit is determined by the detector and the im¬ 
posed upper limit may represent an effort to exclude dijet 
contributions (e.g. p t spectra below 2 GeV/c have been 
interpreted entirely in terms of a thermalized flowing bulk 
medium with no jet contribution [HI)- The lower limit 
has important consequences for ensemble-mean p t and pt 
fluctuation analysis. No jet contribution to p-p or Pb-Pb 
fluctuation data is acknowledged. 

p-p fluctuation data reported at several LHC ener¬ 
gies and corresponding to non-single-diffractive (NSD, 
low-multiplicity) conditions show no significant collision- 
energy dependence over a large energy interval. The p-p 
data are said to exhibit a clear power-law dependence 
on charge multiplicity n c h but deviate from a linear- 
superposition reference oc 1/n^. Several theory Monte 
Carlos (MCs) show “qualitative agreement” with the p-p 
fluctuation data, quite different from the qualitative dis¬ 
agreements for pt{n c h) data reported in Ref. [ 151 . 

Conjectured p-p (p t ) fluctuation mechanisms are reso¬ 
nance decays, jets and quantum correlations. To account 
for those “conventional mechanisms” in Pb-Pb collisions 
p-p results are assumed as a reference, p-p data provide 
a “model-independent baseline....” Nontrivial results in 
A-A would then be signaled by “modification of the fluc¬ 
tuation pattern with respect to the p-p reference.” 

“Fluctuations [of ( p t ) in Pb-Pb] were found to decrease 
with collision centrality, as generally expected in a dilu¬ 
tion scenario caused by superposition of partially inde¬ 
pendent particle-emitting sources.... Deviations from a 
simple superposition scenario have been reported.” The 
linear-superposition reference (for \fC/pt) is again as¬ 
sumed to be a 1 l n °ch- Relative to the p-p power-law 
trend Pb-Pb peripheral data follow a similar power law, 
but more-central Pb-Pb data first rise sharply relative 
to the p-p trend and then fall off more slowly. That is 
noted as a remarkable correspondence given the major 
disagreements for p t data as in Refs. mm- 

The Pb-Pb results are said to be consistent with string 
percolation or the onset of thermalization and collectiv¬ 
ity. No critical behavior is evident. There is possibly evi¬ 
dence for initial-state density fluctuations. Several Monte 
Carlos are in “qualitative agreement” with Pb-Pb data. 
HIJING, a model of A-A collisions, is said to follow a 
linear-superposition power-law reference and is inconsis¬ 
tent with the Pb-Pb fluctuation data. Qualitative agree¬ 
ment is reported between Pb-Pb data and Au-Au data 
and with Monte Carlo models that incorporate collective 
phenomena. It is concluded there is no significant energy 
dependence of Au-Au or Pb-Pb fluctuation data. 
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III. STATISTICAL ANALYSIS METHODS 

Different statistical measures applied to the same par¬ 
ticle data may lead to contradictory physical interpreta¬ 
tions. Does that mean statistical analysis is arbitrary, 
that collision mechanisms cannot be inferred from parti¬ 
cle data? Resolution of ambiguity requires detailed com¬ 
parison of measures and the requirement that a valid de¬ 
scription must confront all analysis results consistently. 

I compare several analysis methods in the context of 
conventional statistics, including the central limit the¬ 
orem and Pearson’s correlation coefficient. To simplify 
algebraic relations I eliminate cumbersome summation 
notation where possible and introduce a compact and 
self-consistent symbol set based on common usage. 

The central limit theorem (CLT) provides a basic ref¬ 
erence for fluctuation analysis. The CLT asserts that 
for certain conditions—uncorrelated samples from a fixed 
parent process—certain moments of the sample popula¬ 
tion are invariant under scale transformations mm- 

A. Basic statistical quantities 

A random variable (RV) represents a set of samples 
from a parent process (e.g. density distribution or se¬ 
quence of physical events such as nuclear collisions). I 
assume some detector angular acceptance (A?;, A</>) on 
pseudorapidity r/ and azimuth (j> that may be partitioned 
into some number of bins M. The basic event-wise RVs 
for the present study are Pt and n c h representing sums 
over those charged particles falling within an angular- 
acceptance bin in a collision event (the entire detector 
acceptance or some fraction thereof). We take n c h —» n 
below to lighten the notation. The ensemble-mean bin 
pair number is then 

_ 1 N evt 

n(n-l) = —— n k {n k - 1) (1) 

evt k =1 

= n 2 + nAcr 2 , 

which defines the normalized number variance difference 
= (ct 2 — n)/n as a per-particle measure of number 
fluctuation excess relative to a statistical reference. For 
number fluctuations cr 2 re f = n is the Poisson reference 
variance. Generally Act 2 = (a 2 — cr 2 re f)/n is a per- 
particle measure of variance excess for RV x (consistent 
with Pearson’s correlation coefficient defined below). A 
compound RV such as ( p t } = (Pt + 6P t )/(n + Sn ) (where 
6x is an event-wise deviation from the ensemble mean) 
is a complex statistic with fluctuations represented by a 
series of variances and covariances of elementary RVs. 

In the data analysis and discussion that follows I as¬ 
sume that multiplicities are integrated within one unit of 
r), so for instance dn c h/drj ss n c h/Ar] -A n c h, but in the 
figures the density ratios are made explicit. Multiplicity 
variables are event-wise random variables n x with means 


n x . To simplify notation I omit the bars on multiplicity 
variables unless there is ambiguity. 


B. Variance-based fluctuation measures 

I assume that fluctuating RV x follows a peaked distri¬ 
bution with mean value x and characteristic r.m.s. width 
cx (linear dispersion measure). Variance ct 2 is a conven¬ 
tional fluctuation measure. The single-particle p t vari¬ 
ance (x —» pt) for a specific acceptance bin is 


< = ( 2 ) 

The variance of bin-sum P t given event-ensemble means 
n and p t is Op = ( P t — npt) 2 . The conditional variance 
of Pt given event-wise bin multiplicity n is 

a p t \n = ( Pt-npt ) 2 (3) 

= (J 2 Pt -2pta 2 p tn +p 2 t (T 2 n , 

where pt = Pt/n, and Op measures the n-Pt covariance 
that prompted the introduction of measure <f> Pt (defined 
below) to study equilibration in A-A collisions [20]. 

The “correlator” C [T3] [T5] defined below is a ratio of 
means. Its numerator B can be re-expressed more simply 
as a variance difference 

n,n— 1 

B = ^2 (pt,i - Pt)(pt,j - Pt) ( 4 ) 

i 

= (P t - npt) 2 - n((p t - pt) 2 ) 

B = ° 2 p t \n-™ 2 p t , 

which is zero for CLT conditions (stationary parent pro¬ 
cess and no significant two-particle correlations) [101 US ■ 

B has been described in terms of covariances [mm 
but the usage is misleading. Whereas the structure of B 
expressed in the first line may suggest a covariance repre¬ 
senting correlations (between what two quantities?) the 
reality is a variance difference (as expressed in the third 
line) representing a conditional fluctuation excess and de¬ 
scribing a single RV. The algebraic relation between fluc¬ 
tuations and correlations was established in Ref. |7j. B is 
a specific example (over a limiting scale interval) of total 
variance difference AE 2 ,^ defined in Refs. 

C. Pearson’s correlation coefficient 

A covariance describes the relation between two dis¬ 
tinct RVs, for instance event-wise sums from different ac¬ 
ceptance bins. In some cases a covariance may be related 
to (normalized by) marginal variances of the individual 
RVs. The prototype is Pearson’s correlation coefficient. 
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Pearson’s normalized covariance or product-moment cor¬ 
relation coefficient for a joint event distribution on two 
random variables x a and x b in separate bins (a, b) is [191 


Tab 


g a6 e 


[- 1 , 1 ]. 


( 5 ) 


the covariance a 2 b normalized by the geometric mean of 
marginal (single-bin) variances as a normalization fac¬ 
tor. Marginal variances a 2 may represent actual marginal 
data projections, a mixed-event data reference or an ide¬ 
alization such as a Poisson reference assuming factoriza¬ 
tion of the joint distribution. The geometric mean in the 
denominator implies factorization assuming a reference 
with no a-b correlations (CLT conditions). 

In the present context a bin-sum P t covariance between 
bins a and b can be defined as 


Bab = (Pt - npt) a (Pt ~ np t )b, ( 6 ) 

and the pair-number reference becomes 


D. T fluctuations and the thermodynamic analogy 

Some proposed fluctuation measures are motivated by 
a priori physical assumptions relating to a temperature 
narrative where it is assumed that the event-wise ratio 
Pt/n may serve as a proxy for temperature T = E/n 
in nuclear collisions by analogy with thermodynamics, 
with temperature fluctuations estimated by crj, ~ a lp t )' 
For central-limit conditions one expects <x'f Pi \ —► a' p Jn. 
providing a reference for detecting “non-statistical” tem¬ 
perature fluctuations via difference a 2 p ^ — cr 2 t /n. In 
the large-n (thermodynamic) limit all fluctuations should 
then decrease toward zero as an apparent manifestation 
of thermal equilibration that is assumed to occur at some 
level, at least within the interval p t £ [0.15,2] GeV/c 
thought to exclude dijet contributions El- 

Event-wise mean p t is a ratio of RVs represented by 

(Pt) = ~ ( 11 ) 

n 


i(n — 1) 


n a n b 


n a n b 


( 7 ) 


consistent with CLT conditions and factorization. Pear¬ 
son’s correlation coefficient for fluctuations of sum p t in 
each of bins a and b given the bin multiplicities is then 


r a b = 


Bab 


\J n a a p t ,a 




( 8 ) 


For single bin a the equivalent of Eq. ® is 


B a a (Pt Tlpt ) 2 77. a (T ’t B a , 


( 9 ) 


the marginal variance difference defined in Eq. Q, and 


B a 




( 10 ) 


is a per-particle ratio measure analogous to Acr^ in 
Eq. Q. Within a limited angular acceptance n a and 
cjp tj0 may be assumed constant across the bin system 
and 


Pt,a 


n(J pj for all bins. 


Pearson’s correlation coefficient serves as a template 
for fluctuation measure design. One must decide in each 
case what is the appropriate reference for a variance or 
covariance depending on the context and the hypothe¬ 
sis to be tested. The per-particle format appearing in 
Eqs. ^ and (10) is a placeholder [as opposed to per- 
pair measures such as Eq. (12) and Eqs. ( [T7| ) and (18) 
below]. Fluctuation variations relative to nucleon partic¬ 
ipant number N part or N-N binary collision number N bln 
may be studied with a simple rescaling. Each of n and 
Pt have their own mean-value and fluctuation systemat- 
ics that can and should be studied first in isolation, not as 
ratios that may obscure underlying physical mechanisms. 


for each event of an event class. Note that the symbol 
(pt) was used to represent the ensemble mean in a previ¬ 
ous LHC analysis [.15] (referred to as “inclusive” (pt)), a 
quantity denoted by p t in the present study. In Ref. [Tj 
the event-wise mean is denoted by MEbE(pt) and the 
event-ensemble mean is denoted by M(p t ) following sym¬ 
bol definitions introduced in Ref. [9j. 

The variance of event-wise mean p t is then 

< 4 ) = (M-Pt) 2 ( 12 ) 

f (Pt~ npt) 2 \ 

l n 2 }' 


The proposed per-pair fluctuation measure <x 2 pt ^ — a p t / n 
is problematic for several reasons: (a) The particle mul¬ 
tiplicity for most collisions is small, with large relative 
fluctuations leading to significant but unrepresented co- 
variances as systematic biases [TO] ; (b) per-pair measures 
include an extra factor 1 /h compared to per-particle mea¬ 
sures consistent with Pearson’s correlation coefficient; 
and (c) collision mechanisms other than “temperature” 
variation may produce p t fluctuations that would be mis¬ 
represented by the variance measure in Eq. (12) BM- 


E. Summary of pt fluctuation measures 

When referring specifically to the event-wise mean (RV 
ratio) and its fluctuations I employ symbol (pt)- More 
generally I refer to p t fluctuations. Several p t fluctuation 
measures have been defined at the SPS and RHIC. 

The $ Pt measure [9] (NA49 collaboration) was moti¬ 
vated by a previously-observed n-Pt covariance in p-p 
collisions as a possible indicator of equilibration in A-A 
collisions and is defined as a difference between r.m.s. 
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quantities J2H] 

$ p t = y/vp^Jn-Vpt ( 13 ) 
(v/^ti"/” + CTpt ) $pt = 

~ 2cr pt $> Pt . 

An initial RHIC pt fluctuation analysis .115] (STAR 
collaboration) defined a related r.m.s. measure Acrp t | n 
(denoted by A u Pt -. n in [10] ) to permit comparisons with 
<f> Pt but adopted a variance-based measure system rather 
than r.m.s. quantities 

Aa P t \u = n ((Pt) ~ Pt) 2 ~ a p t (14) 

= B/n 

= ‘^ ( 7 Pt A( ^Pt\n 

~ 2a Pt . 

The PHENIX collaboration defined a similar r.m.s. mea¬ 
sure [HI 


F Pt = 


a (pt) 


(Tn 


Tl 


A a 


'pt 

Pt\n 


J Vt 


(15) 


where a Pt /y/n approximates a mixed-event reference. 

In subsequent STAR analysis 12 M variance difference 
A°p t | n was redefined as the statistically simpler form 


A( 7 Pt\r 


= B/f 


(16) 


equivalent to a 2 a r a in Eq. (10). Those related per- 


particle measures are based on event-wise bin-sum P t 
fluctuations conditional on bin-sum multiplicity n. De¬ 
nominator n in Eq. (161 serves as a placeholder for quan¬ 


tities such as N par t or Nbi n that better probe fluctuation 
excess in relation to conjectured collision mechanisms. 

In an alternative STAR fluctuation analysis 11 .'i| a (p t ) 
fluctuation measure was based on the assumed tempera¬ 
ture narrative described in Sec. cud] 


2 

a p t ,dyn 


B 

n(n — 1) 


(17) 


’(pt) 




A ° 2 p t \n/( n - 


Aa 2 


covariances, 


which includes a dominant 1/n trend even if the underly¬ 
ing correlations do not change with system size or scale. 

The main subject of the present study is the recent 
ALICE (p t ) fluctuation analysis JTj employing a measure 
similar to the alternative STAR measure in Eq. © 


n(n — 1) n 2 

= Acr P t |n/(™ + Ao£). 


(18) 


C is referred to as a “correlator” intended to represent co- 
variances averaged over particle pairs, but quantity B is 
actually a variance difference as demonstrated in Eq. Q , 
not a covariance. Note that n + Acr^ —> n{\ + A a^Jn) 
and the second term (a per-pair fluctuation measure) be¬ 
comes small compared to 1 for larger multiplicities. 

The measure of (p t ) fluctuation “strength” actually de¬ 
fined in Ref. T] is the r.m.s. quantity 


Vc 7 

p't 


B' 

n ,2 p' t 2 


1/2 


B _^/ 2 

Pt 


(19) 


where primes indicate statistical quantities derived from 
accepted particles only, not corrected to a full p t accep¬ 
tance extending down to zero. Some consequences of 
a low-p t cut limiting the pt acceptance are discussed in 
Sec. [V] and distinctions are maintained between primed 
and unprimed quantities in what follows. 

Ratios of statistical quantities, even mean values, may 
obscure information or combine systematic effects in con¬ 
fusing ways (e.g. variance difference B divided by p t or 
a 2 t ). Fluctuation systematics reflecting strong deviations 
from statistical references may be obscured by such ra¬ 
tios. Systematic trends of primary RVs representing ex¬ 
tensive quantities (e.g. n c h, Pt, B) should be considered 
in isolation before any ratios are introduced. 


IV. COLLISION MODELS 

Basic collision models include the Glauber model of 
A-A collision geometry, the two-component model of p-p 
and A-A hadron production near mid-rapidity and the 
specific TCM for hadron yields vs A-A centrality. 


A. A-A collision-geometry model 

Interpreting statistical trends vs nucleus-nucleus (A-A) 
collision centrality requires a model for the A-A colli¬ 
sion geometry. The Glauber model can be used to relate 
certain A-A geometry parameters to charge multiplicity 
n c h (integrated within some angular acceptance) via the 
A-A total cross section [21]. Glauber-model parameters 
include N part / 2, the number of nucleons (from one nu¬ 
cleus) that participate in at least one N-N encounter, 
and Nbi n , the number of N-N binary encounters (based 
on assumed N-N interaction cross section ctnn)- In re¬ 
lation to the TCM (next subsection) the preferred cen¬ 
trality measure is mean participant-nucleon pathlength 
v = 2 Nbin/N part . Also relevant are estimated impact 
parameter b and A-A transverse-overlap eccentricity e. 

One should note that p-p collisions are not described 
by the eikonal approximation basic to the A-A Glauber 
model wherein Nbi„ oc N p ^ t . For p-p collisions the ana¬ 
log to N par t is n s oc the number of participating low-x 
gluons and the analog to Nb in is oc n 2 s oc the num¬ 
ber of participant gluon-gluon encounters. In either case 
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the number of participant binary collisions oc N) nn or n 2 
predicts the nominal TCM dijet production rate [221 . 

B. Two-component model — TCM 

The p-p TCM was initially based on phenomenological 
spectrum analysis [ 25 ] but has since been related to dijet 
production and QCD theory [ 22 ] [ 29 ] [ 30 ] . The TCM has 
been quite successful in describing a variety of RHIC and 
LHC p-p , p -A and A-A data mmmm- TheTCMfor 
yields, spectra and correlations is based on the assump¬ 
tion that hadron production near mid-rapidity proceeds 
via soft (projectile-nucleon dissociation) and hard (dijet 
production) mechanisms assumed to be linearly indepen¬ 
dent. The soft component is assumed to be universal, 
the same for all systems and collision energies. The hard 
component representing minimum-bias (MB) dijets fol¬ 
lows a non-eikonal trend for p-p collisions and an eikonal 
trend for A-A collisions with larger A [22]. The trend 
for p-A collisions is not known a priori but may involve 
a smooth transition from p-p (N-N) to A-A [II] . 

For produced quantity X (e.g. extensive variable n c h , 
P t or B) the two-component model (TCM) is expressed 
by 


The energy dependence of the TCM is based on the 
empirical relation n s oc log(v / s/10 GeV) for the soft com¬ 
ponent (number of low-* gluons near rj = 0) 22 ] . The 
multiplicity hard component nh (~ dijet fragment yield) 
for p-p collisions then scales as rih oc log 2 (y / s/10 GeV), 
where intercept 10 GeV is inferred from dijet systemat- 
ics oanu- For X = P t in p-p collisions soft component 
p M is observed to be constant over a large energy inter¬ 
val, whereas hard component pt.h is linearly related to the 
MB jet spectrum width oc log(-v/s/3 GeV) [14] (and see 
Fig. i right). Jet-related p t angular correlations inferred 
from per-particle pt fluctuation measure Aa 2 p ^ n = B/n c h 
are also observed to scale oc log(y / s/10 GeV) [5]. 

C. Hadron production model 

Figure [l] (left) shows hadron yields obtained directly 
from identified-hadron spectra (points [ 23 ]) for 200 GeV 
Au-Au collisions. The dash-dotted line shows the conven¬ 
tional 200 GeV Au-Au TCM with fixed x = 0.095 |24| . 
The dashed line is the GLS reference with x = 0.006 x 
2.5 = 0.015 predicted from p-p spectrum data [ 25 ]. The 
solid curve that describes a smooth transition from one 
limiting case to another is defined by 


X = n s x s + n h x h (\fs) = X s + X h for p-p (20) 
= ( N part /2)X S + Nbi n Xh for A-A, 

where x s and Xh(\fs ) in p-p collisions or X s and 
Xh in A-A collisions are fixed quantities in a linear- 
superposition (LS for p-p) or Glauber linear superposi¬ 
tion (GLS for A-A) model but may deviate from a fixed 
TCM reference for real collisions [e.g. variation of Xh(v) 
with v for more-central A-A collisions]. The argument 
of Xh{\fs) indicates that the hard component of X can 
have an energy dependence relating to the underlying 
scattered-parton spectrum whereas soft component x s 
typically does not. GLS represents eikonal linear su¬ 
perposition of participant N-N pairs in A-A collisions 
whereas LS represents non-eikonal linear superposition 
of participant low-* gluon-gluon pairs in p-p collisions. 

Participant number N part in A-A collisions corre¬ 
sponds by hypothesis to n s in p-p collisions, and Nbi n cor¬ 
responds to rift oc n 2 s . Thus, v = 2Nbi„./N part oc N p ^ t for 
A-A collisions corresponds to rih/n s oc n 2 s /n s = n s for p-p 
collisions. Per-participant scaling of X/n s or 2X/N part , 
proportional to n s or v respectively in the p-p and A-A 
systems, is interpreted as an indication of dijet produc¬ 
tion within a linear-superposition scenario. 

The TCM relation between p-p and A-A collisions can 
be expressed by 

« X„ + (y-l)X' h for A-A, (21) 

■L'part 

where X pp = X s + Xh for NSD or MB p-p ss N-N colli¬ 
sions and X' h represents a modified N-N hard component 
for secondary N-N scatters in A-A collisions m- 


— - n c h = n pp [l + x(v)(v — 1.25)] (22) 

^ part 

x(v) = *0 + *i{l + tanh[(i/— ft' 0 )]}/2. 

For Au-Au data n pp « 2.5, *o = 0.015, *i = 0.08, and 
vq = 2.5 represents the sharp transition from Ref. HU. 


■a 



4 
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FIG. 1: Left: Per-participant hadron production mea¬ 

sured by (2 /N P art)dn c h/dri vs v for 200 GeV Au-Au colli¬ 
sions (points) inferred from analysis of identified-hadron spec¬ 
tra [23j . The dash-dotted line is the conventional TCM with 
fixed x = 0.095 [24]. The solid curve is described in the text. 
Right: Hadron production vs v for 2.76 TeV Pb-Pb collisions 
from Ref. [2^ (solid triangles). The dash-dotted line and solid 
curve are the TCM for 200 GeV Au-Au collisions scaled up by 
factors 1.84 (soft component) and 1.84 2 (hard component) re¬ 
flecting soft multiplicity n s scaling as liRy^s/lO GeV) 1221(271 . 

Figure [I] (right) shows 2.76 TeV Pb-Pb hadron pro¬ 
duction data from Ref. [26] (points) compared to a cor¬ 
responding TCM scaled up from the 200 GeV Au-Au 
trend (solid curve). For 2.76 TeV factor 1.84 « 
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ln(2760/10)/ln(200/10) predicts the expected increase in 
n s ,NSD ~ n pp — > 4.6, scaling with low-a; parton partici¬ 
pants as described in Ref. [22]. The same factor is applied 
to x{v) oc n s . The functional form of x(v) at 2.76 TeV is 
very similar to that at 200 GeV with the exception that 
the sharp transition (ST) in jet structure near v = 3 
first reported in Ref. m is shifted down to v ss 2 at the 
higher energy, as first noted in Ref. 123 - 

Equation (22) is used in this study to relate reported 
n c h values from Ref. [d] to fractional cross sections and 
then to Glauber parameters N part / 2, Nf,i n and v = 
ZNbin/Npart according to the methods in Ref. [211 . 


V. p t TCM FOR p-p COLLISIONS 

The TCM for previously-published pt (ensemble-mean) 
data is required to process the p t fluctuation data in this 
study. We summarize TCM results from Ref. [143 that 
describe accurately the required p-p pt data. Thep* TCM 
for Pb-Pb data is summarized in App. [A] 


A. pt TCM definition 


The pt TCM for p-p collisions is based on total pt de¬ 
noted by Pt integrated over all particles within some an¬ 
gular acceptance. If n c h is the mean total charge inte¬ 
grated within the same acceptance ensemble-mean p t = 
Pt/n ch . Just &s_n c h =n s + n h l assume P t = P tjS + P t ,h- 
The TCM for Pt(n s ) is then ?r s p tjS + nhpt,h and for p t is 


p t (n s ,ffs) 


n s pt, s + n h p t ,h(Vs) 
n s + n h 

pt, s + x(n s )p t ,h(Vs) 
1 + x(n s ) 


(23) 


where x(n s ) = an s with a « 0.006 [25]. n s is obtained 
from n ch by n s = (l/2a)[\/l + 4 an ch - 1] since n ch = 
n s + an 2 S . The two P t TCM components can be inferred 
from pt spectrum TCM model functions or from p t data. 

If the pt spectrum is cut off at some small value (e.g. 
Pt,cut « 0.15 GeV/c) 


Pt( n s,Vs ) 


n 'sPt,s + n h pt,h{y/s) 
n' s + n h 

pt,s + x(n s )pt,h(Vs) 
n'Jn s + x(n s ) 


(24) 


assuming no loss from the hard components of n c h and 
P t . We observe that the product Pf. = n' B p' t s « n s p t}S = 
P M is (and Pt fluctuations are) insensitive to a low-p t 
cutoff because only a small fraction of integrated P t . s 
(and none of Pt,h) is involved whereas a substantial 
fraction of n s may be affected. Thus, only the ratios 
p’ ts « Pt lS /n' s and n' s /n s are sensitive to a low-p t cutoff. 

In this study I correct affected quantities for a low-p t 
cutoff at pt, C ut ~ 0.18 GeV/c based on Ref. [2]. From 


reported (and corrected) n c h in one unit of p I obtain 
n s from n c h as defined above. Then nh = n c h — n s and 
n' s = 0.75n s giving the sum n' ch = n' s + hh■ The corrected 
ensemble-mean p t is pt « {n' ch /n c h) Pt given P t ~ P[. 


B. TCM description of pt data 


Figure [2] (left) shows LHC p' t data from p-p collisions 
at 0.9, 2.76 and 7 TeV (upper points) [T5 . The LHC 
particle data were obtained with a nominal pt, C ut = 0.15 
GeV/c, but multiplicity n c h was extrapolated to zero pt- 
Also included are reference pt data from UA1 (open tri¬ 
angles, open circles [28]) and STAR (solid points [25]) 
obtained by model fits to spectra. The curves are defined 
by Eqs. (23) or (24) with parameters from Ref. US¬ 


's; 

o 



FIG. 2: Left: p t vs n c h for several collision energies (rep¬ 
resented as ( pt) in Ref. [Tl]). The upper group of points is 
pt from Ref. m biased by a low-pt acceptance cut near 0.15 
GeV/c. Right: LHC p' t data from the left panel multiplied by 
factor n' ch /n a that removes the bias from the low-pt cut and 
the jet contribution to the denominator of p't- The unive rsal 
soft component pt,s is then subtracted according to Eq. ( |25| ) 
isolating the P t hard component as x{n s )pt,h{\fs) = Pt,hjris- 


Figure [ 2 ] (right) shows the quantity 

x(n s )Pt,h(y/s) 


^p-' t (v^)-p- M 

n. 


(25) 


where the expression on the right follows from Eq. (241, 
pt, s = 0.385 GeV/c is assumed for all energies and 
x(n s ) = a n s /Ar] with a = 0.0055 for A?/ = 1.6. The first 
term of n' ch /n s = n' s /n s + x(n s ) is adjusted so that the 
various data sets have a common intercept at 0. The val¬ 
ues are consistent with an effective p t , C ut ~ 0.18 GeV/c. 


C. pt energy dependence and relation to MB jets 

Figure [3] (left) shows p t data in the form 

- 7 -T (— P'tiVs) -Pt,s] = Pt,h(Vs) (26) 
x(n s ) \n s ) 

for four energies, where pq s has fixed value 0.385 GeV/c 
for all energies. Most of the p t ^ values fall in narrow hor¬ 
izontal bands, but the significant downturn for smaller 














multiplicities is a real feature of the spectrum hard com¬ 
ponent first observed for spectra from 200 GeV p-p col¬ 
lisions with smaller n c h [25]. The solid curve is a 7 TeV 
parametrization used in the present study. 



FIG. 3: Left: Hard components pt,h (represented here as 
{pt)h from Ref. [14]) isolated according to Eq. ( |26|) . Most 
data fall within narrow horizontal bands reflecting point-to- 
point data variation consistent with the p-p TCM. The solid 
curve is a parametrization of 7 TeV data used in the present 
study. Right: p t p mean values from the left panel plotted 
vs parameter A ymax = ln(-y/s/3 GeV) describing variation 
of minimum-bias jet-spectrum widths with p-p collision en¬ 
ergy m- 

Figure (right) shows pt,h{\fs) mean values from the 
left panel (solid points) vs Ay max = ln(-y/s/3 GeV) from 
Ref. [22] where it was demonstrated that jet spectrum 
widths scale with p-p collision energy as Ay max . The 
right panel reveals that p t .h is linearly related to the MB 
jet spectrum width. That trend is consistent with the 
results of Ref. [22] where it was demonstrated that the 
spectrum hard component that determines ptp is pre¬ 
dicted by folding an ensemble of parton fragmentation 
functions with a MB jet spectrum. The hadron spec¬ 
trum hard-component width should then scale linearly 
with the jet spectrum width, and ptp should have the 
linear relation to A y max demonstrated above. For p-p 
collisions the pt vs n c h systematics compel a jet interpre¬ 
tation for the TCM hard component. The pt soft compo¬ 
nent remains consistent with longitudinal-projectile (nu¬ 
cleon) fragmentation independent of collision system or 
energy. Open symbols are predictions for lower energies. 

VI. LHC p-p p t FLUCTUATIONS 

In Ref. [Ij p-p ( pt ) fluctuation data are reported for 
Vs = 0.9, 2.76 and 7 TeV. Plotted as VC’/p't vs n c h the 
data for different energies are not distinguishable (Fig. 1 
of Ref. [T] ). We first focus on the 7 TeV data and later 
estimate the energy dependence of variance difference B. 

A. p-p pt-fluctuation multiplicity dependence 

Figure [I] (left) shows VC'/p't data (points) vs corrected 
rich from Ref. [T] (Fig. 4 left) plotted in a log-log format. 


A “power law” trend oc 1 /V n c h noted in Ref. Q] is also 
shown (solid line). In Ref. [T] p-p data were fitted with 
a power-law model function oc 1 /n b ch with b ss 0.4. The 
dotted curve is 0.244/n°^ 42 . The dash-dotted curve rep¬ 
resents a fit to HIJING data reported in Ref. |Tj. The 
dashed curve is explained below. Among several plot¬ 
ting formats appearing in the present study this format 
is comparatively insensitive to information in the data. 



FIG. 4: Left: Representation of (pt) fluctuation data from 
Ref. [I] Fig. 1 (points) with a measure that exhibits similar 
trends for three LHC energies. The curves are described in 
the text, mpt' stands for uncorrected p' t biased by a low-pt 
acceptance cut. Right: Data from the left panel transformed 
according to the axis labels but employing corrected charge 
density n c h/Ay. The curves are described in the text. 

Figure[I] (right) shows n c hC'/p' 2 « B/n c hPt- The data 
(points) increase monotonically with multiplicity over a 
30% interval. The dotted curve is the power-law model 
in the left panel suitably transformed to obtain a trend 
oc n/.j/ 6 . The dash-dotted line is A 2 = 0.315 2 from the 
left panel. The hatched band represents ( B pp /n c h)/Pt ~ 
0.0145/0.43 2 = 0.078 with the NSD p-p value B vp defined 
in Sec. |VI B| The dashed curve is explained below. This 
format is more sensitive to information in the data but 
does not reveal the underlying production mechanisms. 

The pt TCM summarized in Sec. |V| demonstrated that 
Pt vs n c h systematics (rapidly increasing for p-p colli¬ 
sions) are mainly controlled by dijet production [14] . The 
denominator P 2 in C'/p' 2 ss B/Pf may then obscure the 
mechanisms driving variance difference B in the numer¬ 
ator. To better isolate underlying p t fluctuation mecha¬ 
nisms one should remove the extraneous factors p' 2 and 
n' 2 h from the ALICE measure to obtain variance differ¬ 
ence B as an extensive pt fluctuation measure. 

Figure (left) shows n' 2 h C ~ B vs soft multiplicity 
n s (points) with n s inferred from n c h as described above. 
The transformation from VC'/p't data to C' is based on 
the p' t TCM described in Sec. |Vj including the 7 TeV 
parametrization oipt,h (solid curve) in Fig. [ 3 ] (left). The 
dashed curve LS through the data is explained below. 
The dotted curve is the “power-law” expression from 
Fig. [3] (left, dotted curve) suitably transformed. 

A striking feature of the B trend is the large values 
for higher multiplicities. Both dijet production and the 
Pt-fluctuations hard component scale as n 2 , increasing 
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70-fold in the n s interval from 5.3 (NSD) to 45 [25]. At 
the upper n s limit of the B data the hard-component 
multiplicity is about 25% of the soft component, and 
several jets (with mean energy ss 5 GeV) may appear 
per collision within Ap = 1 [30j . For n s = 45 Pt t h ~ 20 
GeV/c (vs P tiS « 17 GeV/c) and B ss 2.5 (GeV/c) 2 . 
The large B values are then consistent with p-p collision 
dynamics dominated by dijet production. 



pa 

a 

u 


ra 

<1 

J3 

c 



FIG. 5: Left: Pt\n variance difference B derived from the 
data in Fig. [ 4 ] (left) according to the axis labels. Factor n'f h 
biased by a low-pt acceptance cut cancels the same bias in C. 
Right: Ratio B/n s data (points) exhibiting the LS scaling 
(dashed line) expected for dijets. In both panels the dashed 
curve (line) represents the LS trend (participant-gluon-pair 
linear superposition) in the p-p TCM. 


Figure [5] (right) shows the same data (points) in the 
form of ratio B/n s vs n s . The dashed line is B/n s derived 
from a TCM for B similar to the p t results of Ref. m 


B = n s b s + n h bh = B s + B h 


(27) 


0.21 


as in Eq. (20), with b s = 0.0085 (GeV/c) 2 and bh 
(GeV/c) 2 assuming nh = an 2 with a = 0.0055 (25] . 
That relation represents the TCM reference for p t fluc¬ 
tuations corresponding to n c h systematics in Ref. (23] : 
linear superposition of participant low -2 gluon-gluon col¬ 
lisions within p-p collisions. For P t expressed in the same 


TCM format [numerator of Eq. (23)] the coefficients are 
p t s = 0.385 GeV/c and p t , h « 1.9 GeV/c at 7 TeV. The 
dashed curves in previous panels are Eq. ( |27| back trans¬ 
formed according to the various measure definitions and 
using the n' ch and p' t TCM from Sec. |Vj The solid curve 
in Fig. [4] (right) is the same back transformation employ¬ 
ing a constant value p tt h = 1.9 GeV/c rather than the 7 
TeV pt t h parametrization (solid curve) in Fig. [ 3 ] (left). 

As with the spectrum analysis in Refs. [237125] soft 
and hard fluctuation components plotted in an appropri¬ 
ate format can be isolated via the trend vs n c h or n s . 
The p-p soft-component estimate in the present case is 
the hatched band in the right panel. The linear trend of 
B/n s is additional evidence for dijet production in p-p 
collisions: soft n s represents N part for low-a; participant 
gluons within protons (instead of participant nucleons 
within nuclei), and n 2 represents Nb in for dijet produc¬ 
tion from gluon-gluon binary collisions 0. As inferred 


from increase of ensemble mean p t with n c h in Ref. [14] 
the dominant mechanism for variation of pt fluctuations 
in 7 TeV p-p collisions is apparently MB dijet production. 


B. p-p pt-fluctuations energy dependence 


According to Ref. [T] quantity y/C/p t (and therefore 
C/p 2 ) vs n c h exhibits no significant collision-energy de¬ 
pendence for p-p over a substantial energy interval, im¬ 
plying that B/n 2 ch p 2 is also nearly invariant for those 
conditions or B ~ p 2 (for given n c h ) in terms of energy 
dependence. But the energy dependence of pt is well 
known from Refs. 0 [15] and the energy dependence of 
fluctuation measure B can be estimated accordingly. 

Figure [5] (right) includes B/n s trends (dash-dotted 
curves) rescaled from 7 TeV by p 2 (n s , y/s) ratios to other 
energies (2.76, 0.9 and 0.2 TeV). The substantial energy 
dependence is evident especially for larger multiplicities. 
The figure can be compared with Fig. [ 2 ] (right) showing 
the Pt,h/n s energy dependence. That a pt variance mea¬ 
sure might have an energy dependence similar to that for 
p 2 is understandable if a single underlying clijet mecha¬ 
nism is common to the two cases m- 

For the analysis of Pb-Pb data below it is useful to 
determine the soft and hard components of B pp ( B eval¬ 
uated for NSD p-p collisions) for several energies by esti¬ 
mating the terms in Eq. (27). The result for 7 TeV NSD 
p-p collisions with n c h/Ar\ ss 5.3 is B pp = B s + Bh = 
0.045 + 0.032 = 0.077 (GeV/c) 2 , where Bh represents the 
result for in-vacuum p-p (first-hit or primary N-N) colli¬ 
sions as opposed to B’ h for an average of secondary N-N 
collisions within A-A collisions [13] . The per-particle 


ratio is B pp / 


0.0145 (GeV/c) 2 . The result for 


pp 


2.76 TeV p-p collisions with n c h/Ap ss 4.6 is B. 
0.039 + 0.024 = 0.063 (GeV/c) 2 with B pp /n s « o'.0135 
(GeV/c) 2 . The result for 200 GeV p-p collisions with 
nch/Arj « 2.5 is B pp = 0.021 + 0.007 = 0.028 (GeV/c) 2 
with B pp /n s = 0.0112 (GeV/c) 2 . The B h estimates all 
assume bh = 0.21 (GeV/c) 2 as for 7 TeV. 


VII. LHC Pb-Pb p t FLUCTUATIONS 

We next consider 2.76 TeV Pb-Pb data from Fig. 8 of 
Ref. [I]. Those VC'/p' t data have been divided by “power 
law” model function A(n c h/Ar])~ 05 fitted to correspond¬ 
ing HIJING data to obtain the constant A ss 0.315. 
Primes denote uncorrected quantities derived from par¬ 
ticles falling within a restricted p t acceptance. 


A. 2.76 TeV Pb-Pb vs 200 GeV Au-Au 

Figure [6] (left) shows 2.76 TeV p-p (open circles) and 
Pb-Pb (solid dots) data in the form yj ( dn c h/dp) C'/p' t A 
plotted vs n c htAt] representing Fig. 8 of Ref. jT]. While 
there is a rough correspondence in the magnitude of data 
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from the p-p and Pb-Pb collision systems there is a strong 
quantitative disagreement in the trends, just as observed 
for the p t vs n c h trends in Ref. m- The intersection of 
p-p and Pb-Pb trends near n c h/Arj = 25 is misleading. 
Where a single p-p collision in that interval may produce 
one dijet a single A-A collision includes many N-N colli¬ 
sions each with the same small probability (<C 1) of pro¬ 
ducing a dijet, representing very different physical con¬ 
texts. Whereas B/n s oc n s in Fig. I 5 J (right) B/n c h oc rq,{ 3 
in Fig. [?] (left). Interpretation is also hindered due to the 
r.m.s. nature of the preferred ratio measure yC/p t . Be¬ 
low I transform the Pb-Pb data to alternative formats. 



FIG. 6: Left: Representation of 2.76 TeV Pb-Pb data from 
Fig. 8 (solid points) and 7 TeV p-p data from Fig. 1 (open 
circles) of Ref. [I] normalized to A/y/n c h from a fit to HIJING 
data. Right: Transformation of the Pb-Pb data from the left 
panel according to the axis labels with A = 0.315. The result 
is equivalent to A <Jp t \ n from Ref. [TO]. n c h 0 is an estimate of 
the charge multiplicity corresponding to b = 0 central Pb-Pb 
collisions. The curves in both panels are described in the text. 


Figure [6] (right) shows the Pb-Pb data in the left panel 
replotted as {n ch /A??) C/p t « Aerp t | n vs n ch /n ch 0 for 
comparison with Fig. 2 of Ref. [IQj (STAR) ( n c ho is an 
estimate of the multiplicity corresponding to b = 0 cen¬ 
tral collisions, as in Ref. m)- The published data have 
been corrected via n' ch —» n c h and p' t —> pt according 
to the expressions in Sec. V A The required p' t [n c h) val¬ 
ues were obtained from the parametrization of 2.76 TeV 
Pb-Pb p t data introduced in Ref. m and summarized 
in App. [A] The value for p-p collisions (hatched band) is 


obtained from the p-p NSD value B pp /n s from Sec. VIA 


and the NSD p t value from Ref. El - (■ B pp /n s )/pt = 
0.0135/0.42 = 0.032 GeV/c. The GLS reference (dashed 
curve) is explained below. There is good agreement 
with Ref. m 130 GeV data in general shape. The 
higher-energy data exhibit larger fluctuation amplitudes 
( 1=3 2-fold increase) as expected for a dijet mechanism. As 
with plots on N part this plot on n c h compresses more- 
peripheral data into a small interval and de-emplrasizes 
the important GLS scaling at the left. 


Reference [5] (PHENIX collaboration) reported pt. fluc¬ 
tuations measured by F Pt « Acrp t \ n /a Pt defined in 
Eq. fll5| . The centrality trend of F Pt on N part in Fig. 2 
of ReL [8] corresponds well with the centrality trend of 
Ao> t | n on n c h/n c ho in Ref. [XOj and to the LHC fluctua¬ 


tion data in the format of Fig. [6] (right). The maximum 
value <J Pt F Pt « 0.3 x 0.035 ~ 0.01 GeV/c appears to con¬ 
flict with Acrp t | ra « 0.05 GeV/c reported in Ref. [TU]. But 
total variance scales approximately with angular accep¬ 
tance (e.g. Fig. |TT] left), which for STAR is 2 x 27r = 4-7T 
whereas for PHENIX (in that study) it is 0.7 xn = 0.77r. 
The acceptance ratio 5.7 thus accounts reasonably well 
for the STAR/PHENIX fluctuation data ratio 5. 

Reference [5] also presents the effect of varying the ac¬ 
ceptance upper limit p™ ax (equivalent to a running inte¬ 
gral) in its Fig. 3, where the maximum rate of increase 
occurs just above 1 GeV/c at the mode of the spectrum 
hard component reported in Refs. [23, [25j and consistent 
with MB dijets as the principal source of p t fluctuations. 

Figure [7] (left) shows the Pb-Pb data from Fig. [6] 
(left) converted to per-particle variance difference 
{n ch /Aif)C m Aa 2 p ^ n = B/(n ch /Ar]) according to 
Eq. (18) and plotted vs mean participant pathlength v. 
The general trend is strong increase with centrality. The 
dashed curve is explained below. When data are plotted 
on mean participant pathlength v the GLS data trend is 
apparent for more-peripheral data. The Pb-Pb data fol¬ 
low the GLS reference up to v = 2.5, suggesting transpar¬ 
ent Pb-Pb collisions within that interval m- The 2.76 
TeV p-p reference value is B pp /n s « 0.0135 (GeV/c) 2 . 


o 


£ 

< 



FIG. 7: Left: n c hC « B/n c h = A uF (points) transformed 
from the 2.76 TeV Pb-Pb data in Fig. |6](right). Right: Equiv¬ 
alent data for 200 GeV Au-Au collisions obtained from a study 
of pt fluctuation scale dependence in Ref. [T] . 


Figure [7] (right) shows equivalent data for 200 GeV 
Au-Au collisions from Ref. [3] reporting a study of fluc- 
tuat ion sc ale (bin-size) dependence of p t fluctuations (see 
Sec. VIII). Those data correspond to bin size (A 77, A/>) = 
(1.6, 27r) matching the ALICE detector acceptance. The 
general variation with centrality is similar, but values for 
the higher collision energy are substantially larger as ex¬ 
pected for a dijet production mechanism. The Au-Au 
data appear to follow a GLS reference (dashed curve) up 
to at least v « 3. The curvature of the GLS trend at the 
higher collision energy (left) is greater than that at the 
lower energy because of increased dijet production and a 
larger hard-component contribution to n c h in the denom¬ 
inator of A(jp t | n = B/n c h ■ The 200 GeV p-p reference 
value is B pp /n s « 0.0112 (GeV/c) 2 (derived from 7 TeV 
p-p data in Sec. VIB). We now introduce an additional 
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factor 2n c h/N part to obtain (2 /N part )B, the variance dif¬ 
ference per participant-nucleon pair. 

Figure [8] (left) shows 2.76 TeV Pb-Pb p t fluctuation 
data from Fig. [7] (left) converted to the per-participant 
form vs path-length v. For v < 2.5 we observe centrality 
scaling oc v consistent with the GLS trend as expected for 
dijet production within transparent A-A collisions (fol¬ 
lowing Nfnn scaling exactly). Above that point there is 
a 25% increase relative to GLS until v ss 5 above which 
the data show a reduction to 10% below GLS. That panel 
can be compared directly with Fig. I (right) 


o 



FIG. 8: Per-particle data from Fig. [T] multiplied by factor 
2 n c h/Np ar t to obtain per-participant trends for A-A collisions 
equivalent to B/n s in Fig. [H] (right) for p-p collisions. The 
dash-d otte d lines represent TCM “first-hit” trends if A-A B' h 
in Eq. (281 is replaced by in-vacuum p-p Bh- 


The GLS trend for more-peripheral 2.76 TeV Pb-Pb 
data (dashed line) is described by 


N, 


part 


B = B pp + B' h {v- 1.25) GLS trend (28) 


= 0.057 + 0.115(i/-1.25) (GeV/c) 2 , 


with B pp = 0.057 (GeV/c) 2 and B' h « 0.115 (GeV/c) 2 
derived from the data in Fig. [8] (left). From the 7 TeV p-p 
systematics of Sec. VIA B s = 0.039 (GeV/c) 2 and Bh = 
0.024 (GeV/c) 2 . The sum B pp = 0.063 (GeV/c) 2 derived 
from p-p data is consistent with intercept 0.057 (GeV/c) 2 
from Pb-Pb data within the data uncertainties. The p-p 
“first-hit” hard component Bh increases about 4.8-fold to 
B' h for subsequent N-N encounters (“wounded-nucleon” 
effect). Equation (28) is reverse transformed to obtain 
Pb-Pb GLS (dashed) curves in previous figures. 

The hard component B' h for secondary N-N collisions 
in Pb-Pb collisions is much larger than Bh for in-vacuum 
p-p collisions. The dash-dotted line indicates the GLS 
trend predicted from p-p data. Secondary N-N collisions 
within A-A collisions appear to produce dijets at a much 
higher rate (B' h /B h ~ 4.8) than isolated p-p collisions, 
consistent with results from the p t study of Ref. El- 

Figure [S] (right) shows 200 GeV Au-Au p t fluctuation 
data from Fig. [7] (right) converted to per-participant form 
vs mean participant path-length v. For v < 3 centrality 
scaling oc v is consistent with GLS (transparent Au-Au 
collisions). The GLS description for more-peripheral 200 


GeV data is 
2 


N 


B = 0.028 + 0.03(i/-1.25) (GeV/c) 2 , (29) 


part 


with intercept B pp = 0.028 (GeV/c) 2 compared to p-p 
B s = 0.021 and B h = 0.007 (GeV/c) 2 (from Sec.|VIB) 


and B' h « 0.03 (GeV/c) 2 . An accurate B pp value for 
200 GeV Au-Au collisions is thus derived from 7 TeV 
p-p data. Ratio B'JB h = 0.03/0.007 = 4.3 for 200 GeV 
Au-Au compares with 4.8 for 2.76 TeV Pb-Pb. 

Reference [I] reports a “strong reduction of fluctua¬ 
tions” for more-central Pb-Pb collisions, referring to the 
data summarized in Fig. [6] (left) of the present study. 
The decrease is associated with possible “onset of ther- 
malization and collectivity.” However, decreases in other 
plotting formats are much less as in Figs. [7] and [8] The 
present study suggests that most of the decrease for quan¬ 
tity C"/p( 2 is due to the increasing dijet contribution to 
P t 2 in the denominator of equivalent ratio B/Pf. 

While some fraction of the decreases relative to GLS 
above v = 4.5 in Fig. [8] might be attributed to “jet 
quenching” two other explanations are possible: (a) As 
noted, fluctuation amplitudes correspond to integration 
of underlying angular correlations including the SS 2D 
peak attributed to MB dijets m- That peak becomes 
elongated on p in more-central collisions. For the STAR 
and ALICE TPCs an increasing fraction of the SS peak 
extends beyond the detector p acceptance with increasing 
centrality, (b) In each of the more-central A-A collisions 
many dijets appear within the ALICE or STAR TPC ac¬ 
ceptance, resulting in overlap of dijet structure on the 
space ( p , (f>) and failure to fully resolve individual jets - 
a source of inefficiency for fluctuation measurements. 


B. A-A pt-fluctuations energy dependence 


The underlying energy dependence associated with the 
TCM for pt fluctuations from p-p collisions is n s oc 
log(A/s/10 GeV) as reported in Ref. [52]. For example, 
the n s ratio factor is 2.18 for 7 TeV and 1.84 for 2.76 
TeV relative to 200 GeV. The basic n s (oc participant 
low-a; gluons) energy trend permits prediction of NSD p-p 


values for B pp and B pp /n s as summarized in Sec. VIB 


Given the basic logarithmic n s (y/s) dependence the nom¬ 
inal TCM relation between p-p and A-A fluctuations is 

B pp {Vs) = n s b s + an 2 s b h (y/s) = B s + B h for p-p (30) 
B(v, y/s) = B pp {\fs) + B' h (v, \fs){v — 1.25) for A-A, 

with n s in the first line given the NSD value and a ss 
0.006. For proper comparisons the observed n c h must be 
relat ed acc ura tely t o inferred n s and v as described in 
Secs. IV C and V A From Ref. [21] the predicted v value 
corresponding to NSD p-p or N-N averaged over an A-A 
collision is 1.25 as in Eq. (30). In that context I consider 


the energy dependence of p t fluctuations in A-A collisions 
where there are two issues: (a) comparison of B pp from 
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p-p collisions to an equivalent value for peripheral A-A 
collisions and (b) the energy dependence of B' h for A-A 
collisions in contrast to Bh for in-vacuum p-p collisions. 

Regarding point (a), results in Figs. [7] and [s] indicate 
that R pp ( v / s) values (hatched bands) predicted for two 
lower energies from 7 TeV p-p data agree with the corre¬ 
sponding B(u, y'SAriv) (evaluated at v = 1.25) within the 
data uncertainties. The comparison was made without 
adjustment of fluctuation data or A-A centrality mea¬ 
sures. The pt fluctuation data were obtained by two col¬ 
laborations with different detectors and methods. The 
p-p extrapolation was made assuming bh(\/s) is indepen¬ 
dent of energy, although there is reason to expect a mono¬ 
tonic relation similar to that for p t ,h from Ref. [2]. The 
Bh estimates for the lower energies should then be too 
large, but the relative effect on the B pp estimate is minor 
because of the quadratic decrease of Bh with NSD n s . 

Regarding point (b), in Fig. [8] 13/ = 0.115 and 0.03 
respectively for 2.76 TeV and 200 GeV with ratio 3.8. 
The corresponding k^v^s) ratio is (1.84) 2 « 3.4. If 
B' h oc Bh for both energies we would expect the bh ratio 
to be 3.8/3.4 ss 1.1, consistent with only slight variation 
of bh(y/s) with energy where I have assumed none. That 
result is also consistent with approximately the same ra¬ 
tio B'JBh ~ 4.5 at both 2.76 TeV and 200 GeV. The 
ratio itself remains unexplained. 

The collision-energy systematics for p t fluctuations 
thus strongly suggest that almost all hadron production 
arises from participant low-a; gluons following a simple 
QCD logarithmic energy trend for n s , either directly as 
the soft component or via large-angle scattering to dijets 
as the hard component. Yields, spectra, fluctuations and 
correlations are described quantitatively by the TCM. 
Whereas pt fluctuation data from p-p collisions and more- 
peripheral A-A collisions follow a LS or GLS TCM ref¬ 
erence accurately, the jet-related data from more-central 
A-A collisions deviate from the GLS reference quantita¬ 
tively and may provide insight on A-A jet modifications. 


C. Pb-Pb pt fluctuations vs Monte Carlos 

In Ref. pQ two theory Monte Carlos are compared with 
the Pb-Pb data. HIJING j34j is an A-A model based 
on PYTHIA ;35] (initial-state strings and minijets) plus 
Glauber linear superposition. AMPT PH is a transport 
Monte Carlo based on HIJING initial conditions plus par- 
ton and hadron rescattering in the default version. In the 
string-melting version initial-state partons are combined 
via coalescence to form hadrons. 

Figure [9] (left) repeats Pb-Pb data and results for the 
two Monte Carlos from Fig. 8 of Ref. [IL It is noted 
that HIJING data in the form J(n ch /Ar/)C'/pf t is es¬ 
sentially constant with A ss 0.3 except for the most- 
peripheral points. Default AMPT increases much more 
rapidly than the Pb-Pb data, whereas “string melting” 
AMPT increases much more slowly than the data. The 
HIJING results are discussed further below. 



FIG. 9: Left: Representation of Pb-Pb data (solid points) 
from Fig. 8 of Ref. [T[ including default AMPT (open circles) 
and HIJING (open triangles) Monte Carlo results. Right: 
The Pb-Pb data and AMPT MC results transformed to per- 
participant format (2 /N par t)B. 


Figure [9] (right) shows the same data in the form 
(2 /N part )B measuring the Pt\n variance difference per 
participant pair. Default AMPT rises to twice the Pb-Pb 
data for more-central collisions. Since AMPT relies on 
HIJING for its initial conditions, and HIJING shows 
negligible centrality dependence (see discussion below), 
one can ask what rescattering process can produce such 
large pt fluctuations. The dashed curve represents Nhi n 
scaling of dijet production, the principal QCD mecha¬ 
nism for transporting longitudinal projectile momentum 
to transverse phase space in elementary collisions. That 
the Pb-Pb data exceed that level by 25% in more-central 
collisions is already notable. For reasons given below HI¬ 
JING transport is much less than required by the data. It 
is not clear how AMPT compensates for that deficiency. 

With “string melting” enabled the AMPT fluctuation 
amplitude is much smaller than the Pb-Pb data, presum¬ 
ably for the same HIJING input. What happens to the 
transverse momentum manifested by the default version? 
Since HIJING was formulated to model minijet produc¬ 
tion as described quantitatively by QCD theory why isn’t 
dijet production a basis for discussion of results? 

Figure[l0|(left) shows 2.76 TeV Pb-Pb fluctuation data 
intheformH/(n c / l /Ar 7 ) = A(Tp^ n from Fig. jTj (left). The 
HIJING trend (open circles) is obtained from the data 
summary in Fig. [ 9 ] (left, open triangles) by first squaring 
those data then multiplying by A 2 = 0.315 2 and fixed 
p' 2 = 0.5 2 (GeV/c) 2 assuming negligible p t variation. 

Figure 10 (right) shows 200 GeV Au-Au data (solid 
points) from Fig. [ 7 ] (right) reported in Ref. [4]. The open 
points are from a scale (bin-size) analysis of pt fluctua¬ 
tions (see Sec. [Vm| from HIJING for 200 GeV Au-Au 
collisions [6] corresponding to the ALICE detector ac¬ 
ceptance (1.6,27r). The quantitative relation to the 2.76 
TeV HIJING prediction is notable: HIJING B/n c h en¬ 
ergy variation follows the data B pp /n s trend. 


In an earlier study it was determined that p t fluc¬ 
tuations from HIJING are dominated by dijet produc¬ 
tion 331. Disabling dijets in HIJING produced a five-fold 
reduction in Aa Pt \ n for 0-5% central Au-Au collisions 
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FIG. 10: Data and trends from Fig. [7] compared to HIJING 
MC results. The MC results in the left panel are transformed 
from Fig. [ 9 ] (left) assuming a fixed value p' t = 0.5 GeV/c. The 
MC results in the right panel are obtained from Ref. [6]. 


equivalent to an eighteen-fold reduction in B. Related 
angular-correlation studies show that HIJING correlation 
structure above 0.5 GeV/c is negligible with jets disabled. 
We conclude that almost allp t fluctuations from HIJING 
arise from (mini)jets. 

An explanation for the HIJING fluctuation centrality 
trend is provided in Sec. VIII I of Ref. El relating to 
per-particle amplitude variation of the jet-related SS 2D 
peak in 200 GeV Au-Au number (as opposed to p t ) angu¬ 
lar correlations. The HIJING centrality trend is nearly 
constant while the Au-Au data trend shows an eight¬ 
fold increase from the p-p value. The difference arises 
from two sources. The HIJING MC produces a hard- 
component (jet-related, binary-collision scaling) hadron 
yield per binary collision that is seven times the p-p yield 
and 1.6 times the more-central Au-Au yield, but the num¬ 
ber of jet-correlated hadron pairs is 60% of the p-p equiv¬ 
alent and only 20% of the more-central Au-Au equivalent. 
Thus, a per-particle correlation measure such as At 7 p t | n 
(number of correlated pairs / number of particles) may 
increase rapidly for Au-Au data but not at all for HI¬ 
JING data. The same trends apparently determine pt 


fluctuations for 200 GeV and 2.76 TeV as in Fig. 10 


The underlying difference is the PYTHIA model [32] 
for dijet production in p-p collisions that assumes an 
eikonal approximation inconsistent with p-p spectrum 
data m and a default jet spectrum lower limit po = 2 
GeV for HIJING compared to the observed 3 GeV [22] . 
resulting in a large excess of very-low-energy partons. 


VIII. p t FLUCTUATIONS vs CORRELATIONS 


A. Total-variance scaling and fluctuation inversion 

An early p t fluctuation study of RHIC data reported in 
Ref. [TO] related the concept of total variance introduced 
in Ref. [15] to analysis of fluctuation scale dependence. 
The ensemble-mean total variance for event-wise pt and 
particle-number distributions on space x with M bins of 
scale (bin width) Sx within acceptance Ax is 

M 

Zp t \ n {Ax,5x) = ^2[P t (Sx) - n(Sx)p t ]l, (31) 

a—1 

where P t (Jx) and n(<5x) are bin sums and n c h is the mul¬ 
tiplicity in acceptance Ax. The CLT is equivalent to the 
statement that for certain conditions (independent sam¬ 
ples from a fixed parent process) the total variance is 
scale invariant. In general the total variance approaches 
the limit n c h<j^ t at the “single-particle” scale [Sx <C Ax, 
one particle per bin in M(Ax) = n c h occupied bins]. 
The total-variance difference AEp t | n (Jxi, 8 x 2 ) over some 
scale interval [Jxi, 8 x 2 ] is nonzero if CLT conditions are 
not met: The parent process varies from event to event 
and/or the samples (e.g. particle momenta) are corre¬ 
lated. Quantity B defined in Eq. Q is the P t \n total- 
variance difference evaluated over the maximum accessi¬ 
ble scale interval - between the detector-acceptance scale 
and the single-particle scale. 

Figure [Tl] (left) shows the 2 D scale dependence of 
per-particle variance difference (relative to the single¬ 
particle scale) Ao-p t | n (J? 7 , 8 (f)) = AEp t | n (< 5 ? 7 , 8 (j))/n within 
the STAR TPC acceptance for 200 GeV Au-Au colli¬ 
sions .4], The ALICE Pb-Pb fluctuation data reported 
in Ref. |T] correspond to a single point at ( 1 . 6 , 27t) on a 
similar surface representing 2.76 TeV Pb-Pb collisions. 
The scale variation of AEp^/Jr?, 8(j)) has been expressed 
as the running integral of a 2D angular autocorrelation 
on difference variables (t 7 a, 0 a) in the form of an inte¬ 
gral equation. The underlying angular correlations can 
be inferred by inverting that integral equation [7]. 
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The LHC pt fluctuation measurements correspond to 
a single 2D ( 77 , (f> ) bin size (scale) - the TPC detector 
acceptance. In previous studies the scale variation of 
Pt fluctuations was measured for Au-Au collisions and 
directly related to underlying angular correlations mill 
[7'. To establish a physical interpretation for event-wise 
p t fluctuations at the LHC it is useful to review results 
from published fluctuation scaling studies at the RHIC. 


FIG. 11: (Color online) Left: A<r 2 t;rl (GeV/c) 2 distributions 
on scale (8r), 8if>) for 45-55% central 200 GeV Au-Au collisions. 
Right: Corresponding autocorrelations on difference variables 
(f?A, </>a) inferred from data at left by fluctuation inversion. 

Figure [lT] (right) shows the inferred p t angular cor¬ 
relations gjf The result has three basic elements: (a) 
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a same-side (SS) 2D peak, (b) an away-side (AS) ID 
peak on azimuth and (c) a non-jet cylindrical quadrupole 
[cos(2^a) dependence]. The general combination is 
the same as that observed for number angular correla¬ 
tions nn ug, but there are quantitative differences in 
the SS 2D peak structure. Elements (a) and (b) have 
been identified with dijet production in a number of 
ways [H] [531 HH 13011311133] ■ Element (c) might be re¬ 
lated to elliptic flow if that were relevant to nuclear col¬ 
lisions [361439] . 


B. pt correlations vs pt fluctuations 


Figure 12 (left) shows the best-fit SS 2D peak am¬ 
plitude A 2 d (solid points) for 200 GeV Au-Au data vs 
path length v [3]. The peak amplitude increases with 
centrality to a maximum value and then decreases for 
the most central collisions. The monotonic increase for 
more-peripheral collisions, approximately proportional to 
v (dashed curve), is consistent with the binary-collision 
scaling expected for dijet production. A 2 d is closely cor¬ 
related with the amplitude of the AS-dipole component 
of the 2D fit model, as expected for dijet correlations. 
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FIG. 12: Left: Jet-related same-side 2D peak amplitude 
(solid dots) vs path length v from 2D fits to pt angular corre¬ 
lations from 200 GeV Au-Au collisions as in Fig. 11 (right). 
The dash-dotted and dotted lines represent corresponding 
HIJING quench-off and quench-on results respectively from 
Ref. [6|. There is no soft component to those jet-related data. 
Right: Per-particle pt fluctuation amplitude B/n c h as shown 
in Fig. [7] (right) for comparison with correlations, including 
the p-p estimate B pp /n s ~ 0.0112 (GeV/c) 2 from Sec. 


VIB 


Figure [l2] (right) shows pt fluctuations from 200 GeV 
Au-Au collisions as the integral (up to a specific bin size 
or scale) of p t angular correlations evaluated at the AL¬ 
ICE TPC angular acceptance 0). The two panels are di¬ 
rectly related since angular correlations in the left panel 
are obtained by inversion of the scale dependence of fluc¬ 
tuations as represented in the right panel. 

Note that for 200 GeV p t fluctuations and correla¬ 
tions there is no sign of the “sharp transition” appear¬ 
ing near v = 3 in number angular correlations reported 
in Ref. m The substantial increase in jet-correlated 
hadron pairs above the transition point was attributed to 
strong modification of parton fragmentation in the A-A 


environment [20] , The same study indicated that modi¬ 
fied jets still retain almost all of the parton energy and 
hence most jet-related p t , consistent with Fig. |j~2} 


C. Au-Au p t angular correlations 

p t angular correlations can be obtained either by inver¬ 
sion of the scale dependence of pt fluctuations HE] or 
by direct pair counting. The same jet-related correlation 
structures are observed, with minor quantitative differ¬ 
ences in the SS 2D peak structure. Those results indi¬ 
cate that the (pt) fluctuations expected to reveal critical 
fluctuations of temperature near a QCD phase boundary 
are actually dominated by a MB jet (minijet) contribu¬ 
tion h mi nn m m- 

Figure [l~3] (upper panels) shows pt angular correlations 
for (a) 85-95% and (b) 10-20% central 200 GeV Au-Au 
collisions obtained by inversion of p t fluctuation scale de¬ 
pendence [7. Fitted AS dipole and nonjet quadrupole 
components have been subtracted to isolate the SS 2D 
peak structure. Similar analysis of HIJING Monte Carlo 
data supports a jet interpretation for the SS peak [Bj. 




FIG. 13: (Color online) Upper: pt angular correlations for 
(a) 85-95% and (b) 10-20% central 200 GeV Au-Au collisions 
inferred by inverting p t fluctuation scale dependence [4] . AS 
dipole and nonjet quadrupole components of 2D model fits to 
the data are subtracted. Lower: Results for the same colli¬ 
sion systems but pt correlations are obtained by direct pair 
counting rather than fluctuation inversion. Improved angular 
resolution and unfiltered statistical fluctuations are evident. 
Data for all panels include an additional acceptance factor 47 t. 

Figure [13] (lower panels) shows p t angular correlations 
for the same collision systems obtained by direct pair 
counting, confirming the results in the upper panels ob¬ 
tained by fluctuation scale inversion. The SS 2D peak for 
Pt correlations is narrower than that observed for num¬ 
ber correlations. The difference is expected for jet cor- 
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relations, since fewer fragments with larger momenta are 
found closer to the jet thrust axis and more fragments 
with smaller momenta appear at larger angles. 

From these comparisons between Pb-Pb and Au-Au 
Pt fluctuations and between Au-Au pt. fluctuations and 
jet-related pt angular correlations we may conclude that 
there are strong similarities between the 2.76 TeV Pb-Pb 
p t fluctuation data and 200 GeV Au-Au p t and number 
angular-correlation data identified with dijet production. 
Per-participant measures typically exhibit nearly linear 
increase oc v (GLS) for more-peripheral collisions possi¬ 
bly followed by an interval of significant increase above 
GLS (much larger for number than for p t correlations) 
and significant decrease for most-central collisions. The 
GLS trend for p t fluctuations is consistent with dijet 
production (binary-collision scaling) and with jet-related 
structure in both number and pt angular correlations. 


IX. DISCUSSION 


A. LHC fluctuation analysis and interpretations 


The analysis in Ref. [I] is similar to several previous 
studies at the SPS and RHIC motivated by a search for 
critical fluctuations associated with the structure of the 
QCD phase boundary. It is conjectured that event-wise 
mean p t may represent the temperature of a thermody¬ 
namic state formed in A-A collisions, in which case criti¬ 
cal temperature fluctuations near a QCD phase boundary 
or critical point may be reflected in mean-p t fluctuations. 


Given the definition of quantity C in Eq. (18) as a mea¬ 


sure of “nonstatistical” (p t ) fluctuations (excess variance 
relative to a reference) and the context of the temper¬ 
ature narrative \/C ~ ST represents an r.m.s. measure 
of excess temperature fluctuations and pt ~ To- The 
adopted fluctuation measure yfC/pt then emulates ST/Tq 
as an r.m.s. measure of relative temperature fluctuations. 

It is concluded that there is no significant energy de¬ 
pendence of (pt) fluctuations for p-p or A-A collisions 
over a large energy interval. The p-p data appear to 
decrease according to a power-law trend on n c h■ The 
A-A data decrease with centrality (and therefore n c h ) but 
with substantial deviations from the p-p power-law trend. 
The decrease with A-A centrality is said to be expected 
for a “dilution scenario caused by [linear] superposition 
of partially independent particle-emitting sources.” The 
linear-superposition reference is assumed to be 1 /y/n c h- 
Deviations of Pb-Pb fluctuations from the LS reference 
are said to be “consistent with” (a) string percolation or 
(b) onset of thermalization and collectivity. No “critical 
behavior” is observed (but the concept is not defined). 
Initial-state density fluctuations may also contribute. 

The analysis and interpretation of Ref. [Ti rely on a 
number of assumptions that may be questioned. It is as¬ 
sumed that a thermodynamic state with a well-defined 
temperature may be formed in high-energy nucleus- 


nucleus collisions and (pt) is a temperature estimator. 
But when isn’t a thermodynamic state formed: periph¬ 
eral collisions, p-p collisions? If a phase transition were 
relevant to high-energy nuclear collisions then some re¬ 
lated fluctuation excess might arise. But the assumption 
that “nonstatistical fluctuations” in some statistical mea¬ 
sure may necessarily reveal a phase transition is not jus¬ 
tified. Linear superposition of independent subsystems 
is said to result in reduction (dilution) of relative tem¬ 
perature fluctuations. But if most subsystems down to 
individual N-N collisions are independent how does the 
composite system qualify as a thermodynamic state? 

Use of p-p ( p t ) fluctuations as a baseline or reference 
for A-A collisions is proposed but not implemented. Pos¬ 
sible “conventional mechanisms” for (p t ) fluctuations in 
p-p collisions are acknowledged, including jets, but are 
not pursued further. It is often assumed that jets do 
not contribute to pt < 2 GeV/c in A-A collisions [17] 
(possibly what motivates the pt < 2 GeV/c cut for the 
LHC analysis), but there is plentiful evidence that most 
jet-related hadrons appear below 2 GeV/c and within re¬ 
solved jets [30, HD]. Analysis of LHC p t systematics [14] 
shows that variation of p t with n c h is dominated by dijet 
production, contrary to some assumptions of Ref. |pQ. 

Previous pt fluctuation analyses at the RHIC are cited 
(in Ref. t I] as [15-20]) but only [18] appears again in 
the text - motivating measure C in Eq. (2). Angular- 
correlation and energy-dependence results from [19,20] 
revealing a dijet contribution dominating p t fluctuations 
and angular correlations are not mentioned. Although 
Refs. [18] and [20] disagree diametrically about p t fluc¬ 
tuation energy dependence that is not acknowledged. 

Certain Monte Carlo models are said to be in “qual¬ 
itative agreement” with the fluctuation data, but that 
implies quantitative dis agreement to an extent that may 
falsify the models. HIJING is assumed to repres ent lin - 
ear superposition of N-N collisions (but see Sec. IXCl, 
and any difference from HIJING data should then indi¬ 
cate mechanisms unique to A-A collisions. But HIJING 
based on PYTHIA includes an incorrect model for p-p 
collisions, as noted in Refs. mm, and deviates dramat¬ 
ically from an A-A GLS reference based on measurements 
that represents superposition of real N-N collisions. 


B. Fluctuation measure definitions 

Several pt fluctuation measures applied to data at the 
SPS, RHIC and LHC are reviewed in Sec. |III E] The mea¬ 
sure definitions appear to follow two opposing principles: 
(a) a temperature narrative motivating model-dependent 
“ratios of ratios” (intensive quantities) in which cancel¬ 
lations may conceal significant systematic trends and (b) 
a model-independent approach based on extensive quan¬ 
tities (e.g. n c h, Pt) and their statistical properties. The 
primary event-wise RVs are total Pt and total n c h within 
a single fixed angular acceptance or within each of sev¬ 
eral angular bins covering that acceptance. Deliberate 
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choices must be made regarding combinations of those 
RVs to form secondary fluctuation measures. 

Reference pQ applies the term “fluctuation strength” to 
y/C'/pt, but the term could be applied as well to several 
other statistical measures that incorporate n c h and P t as 
RVs including variance difference B. Both C and pt are 
ratios of mean values then combined to form a secondary 
ratio. Large but similar fluctuations in the primary RVs 
may nearly cancel in the square root of a ratio of ratios, 
possibly obscuring significant collision mechanisms. 

One should first measure fluctuation trends for each 
primary RV separately (variances) and in combinations 
(covariances) relative to well-defined statistical references 
(representing CLT conditions). The choice of fluctuation 
measure should be compatible with correlation measures 
given the direct algebraic relation [7|. Data presentation 
involves both y-axis (fluctuation measures) and x-axis 
(system “size” measures n c h, N part , v) choices. A good 
plotting format may test a significant hypothesis (e.g. 
Fig. I right, Fig. [ 8 ]) while a poor format may obscure an 
important data trend (e.g. Fig. [4] left, Fig. [ 6 ] right). 


C. Linear-superposition references 

Reference JQ defines an “independent superposition” 
(of unspecified particle sources) reference as yfC/pt oc 
1 /y/n c h or n c hC/pt oc constant. But n c hC « Acrp t | ra and 
the ratio Aa 2 p ^ n /p 2 would remain independent of A-A 
centrality only if the numerator and denominator happen 
to vary in proportion or each remains constant. The 
latter would hold in the absence of dijet production. The 
quantities would then include only soft components with 
A(T Pt\n B s! n s and pt pt., s - The defined reference 
is thus equivalent to claiming that dijet production does 
not contribute significantly to (pt) fluctuations. 

Linear (independent) superposition of particle sources 
is the basic assumption of LS references for the TCM 
in which (at least) two specific hadron sources are as¬ 
sumed to contribute: The soft component represents 
projectile-nucleon dissociation (longitudinal fragmenta¬ 
tion) and the hard component represents transverse frag¬ 
mentation of large-angle-scattered parton pairs to dijets. 

For a LS reference jet-related contributions to yields 
and spectra [23, 3Sj as well as to fluctuations and cor¬ 
relations should scale with p-p multiplicity as n 2 assum¬ 
ing linear superposition of gluon-gluon binary encounters 
(e.g. Fig. i left) and with A-A centrality as Nun assum¬ 
ing a GLS reference for transparent A-A collisions. Per- 
participant measures should then scale linearly with n s 
for p-p collisions (e.g. Fig. [5j right) and with v for A-A 
collisions (e.g. Fig. [ 8 ]). A MB dijet (average over the MB 
dijet spectrum at given energy) is observed to contribute 
a certain fixed amount to pt to pt fluctuations B 
and to number and p t angular correlations on ( 17 , <j>) [ 3 ]. 

Because each of n c h, Pt , B and jet-related spectrum 
and correlation components has a unique TCM represen¬ 


tation (but with similar forms) an LS reference for any 
one quantity cannot be simply expressed in terms of an¬ 
other, as proposed in Ref. pQ. It is only in the limit of 
no dijet contribution that the surviving soft components 
would be simply related, with n c hC oc constant and pt 
independent of p-p multiplicity or A-A centrality. HI- 
JING with jet production enabled and no jet quenching 
follows that trend approximately, but detailed study of 
HIJING yields and correlations [IT] reveals that HIJING 
overpredicts the hard-component multiplicity while un¬ 
derpredicting hard-component jet correlations (relative 
to A-A data) such that the ratio B/n c h ~ n c hC is nearly 
independent of A-A centrality as in Fig. [To] HIJING is 
thus not representative of GLS scaling in A-A collisions. 


D. The dominant role of minimum-bias dijets 

Reference |lj does acknowledge jets as a possible “con¬ 
ventional” mechanism for particle production. It is im¬ 
plied that jet contributions to p-p collisions might then 
be used to identify equivalent structure in A-A collisions, 
but the only comparison of p-p and Pb-Pb structure in 
that study is Fig. 5 (or equivalently Fig. 8 ) where no jet 
contribution is identified. 

In the TCM context a strong dijet contribution to LHC 
Pt fluctuation data is easily identified for both p-p and 
Pb-Pb collisions. Figure [5] shows variance difference B 
varying with n s in exact accord with the p-p TCM LS 
reference within data uncertainties. The hard compo¬ 
nent (dijets) is observed to dominate p-p pt fluctuations, 
varying in proportion to n 2 representing gluon-gluon bi¬ 
nary collisions. In Fig. [ 8 ] (left) data for more-peripheral 
Pb-Pb collisions again show agreement with GLS scaling 
expected for dijet production in A-A collisions. With in¬ 
creasing A-A centrality measure B exceeds Nf nn scaling 
by about 25%, consistent with persistence of the dijet 
mechanism but with some quantitative modification. 

An analysis of LHC pt vs n c h data for several energies 
and collision systems indicates that p t data trends are all 
accurately described by the TCM. The p-p data follow 
LS scaling across a ten-fold increase of n s . The dijet pro¬ 
duction rate then increases 100 -fold, implying multiple 
MB dijets per p-p collision on average [H] and consistent 
with LHC p-p number angular correlations mm- 

200 GeV Au-Au p t fluctuation data show trends very 
similar to the LHC data represented in this study 00 . 
Those results are in turn consistent with number correla¬ 
tion measurements indicating that dijets are a dominant 
particle production mechanism El and suggesting that 
jet manifestations are very similar at the RHIC and LHC 
modulo a basic QCD log(-^I/10 GeV) scale factor for n s . 

Previous to the present study a wealth of evidence 
for dijet dominance of high-energy nuclear collisions has 
been presented. The MB dijet-based TCM context is 
internally consistent and has been employed to predict 
and explain many experimental results from the RHIC 
and LHC, including (a) systematics of LHC ensemble- 
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mean p t [2] , (b) pt angular correlations HH] , (c) num¬ 
ber angular correlations mmm 135] . (d) trigger- 
associated transverse-rapidity correlations [SHU], (e) Pt 
spectra [231 EH] and (f) jet-related systematics of hadron 
yields mi esi mi. Those results are all in accord with 
measured dijet properties [22] 00] and QCD theory 29j. 


X. SUMMARY 

A measurement of fluctuations in event-wise mean 
transverse momentum denoted by ( p t ) from p-p and 
Pb-Pb collisions at the large hadron collider (LHC) has 
been reported recently. The fluctuation measure denoted 
by VC/pt is motivated by a temperature narrative in 
which collisions attain some degree of thermalization and 
are characterized by a temperature estimated by ( pt ) as 
one property of a thermodynamic state. Excess ( p t ) fluc¬ 
tuations compared to a reference might indicate the pres¬ 
ence of a phase boundary between a conjectured quark- 
gluon plasma (QGP) phase and a hadron-fluid phase. 

It is inferred from the LHC data that ( p t ) fluctuation 
“strength” is nearly independent of collision energy over 
a broad interval for both collision systems. For p-p col¬ 
lisions fluctuations are said to decrease with increasing 
particle multiplicity n c h approximately as a power law 
n ch A - F° r Pb-Pb collisions ( p t ) fluctuations also decrease 
overall, but relative to the p-p trend they increase for 
mid-central collisions and then decrease for most-central 
collisions. The Pb-Pb results are said to be consistent 
with models that incorporate collective phenomena. 

The choice of fluctuation measure \/C/pt from among 
several candidates is a critical step in such data analysis. 
In the present study I review several measures applied 
previously to fluctuation analysis of nuclear collision data 
and describe their algebraic relationships. I identify to¬ 
tal multiplicity n c h and total transverse momentum P t 
(falling within some detector angular acceptance) as the 
basic extensive random variables for the data system, 
with ( p t ) = Pt/n c h as a derived intensive ratio. I in¬ 
troduce variance difference B for P t conditional on n c h 
as a physical-model-independent p t fluctuation measure. 

Given the relation C « B/n 2 h , data from the LHC ( p t ) 
analysis can be converted to other formats and compared 
directly with previous analysis at the relativistic heavy 
ion collider (RHIC). RHIC fluctuation analyses employ¬ 
ing a per-particle variance-difference measure in the form 
B/n c h revealed pt fluctuations increasing strongly with 
Au-Au collision centrality and with collision energy, very 
different from the reported LHC trends. The scale (angle- 
bin-size) dependence of B/n c h was also measured and 
inverted via a standard numerical method to reveal the 
underlying p t angular correlations. Principal features of 
the inferred p t correlation structure were identified with 
minimum-bias jets and were subsequently confirmed by 
correlation analysis based on direct pair counting. 

The LHC p-p (p t ) fluctuation data, when converted to 
measure B , are described accurately by a two-component 


(soft+hard) model (TCM) in which the hard component 
represents minimum-bias (MB) dijets. The p-p TCM has 
been successful in describing yield, spectrum and correla¬ 
tion data at the RHIC and, most recently, ensemble-mean 
pt vs n c h trends from the LHC. The TCM description of 
B vs n c h for LHC p-p collisions confirms that MB dijets 
dominate p t fluctuations for larger event multiplicities. 

The TCM can also be applied to A-A collision data, 
with Glauber linear superposition (GLS) of nucleon- 
nucleon (N-N) collisions within A-A collisions as a ref¬ 
erence. The TCM description of Pb-Pb ( p t ) fluctuation 
data converted to B indicates that p t fluctuations follow 
a GLS trend with collision centrality for more-peripheral 
Pb-Pb collisions (indicating transparency) but deviate 
quantitatively from that trend for more-central collisions. 

The energy dependence of p t fluctuations from p-p col¬ 
lisions measured by B are predicted over a range from 
RHIC to LHC energies by a simple log(s/so) trend con¬ 
sistent with QCD field theory and with measured system¬ 
atics of MB jet spectra. B values from peripheral Pb-Pb 
(« N-N) collisions are quantitatively consistent with the 
values from non-single-diffractive p-p collisions, including 
the energy dependence from RHIC to LHC. In that com¬ 
parison B pp values for RHIC 200 GeV p-p collisions were 
successfully inferred from LHC 7 TeV p-p y/C/p t data. 

The overarching message from the LHC p t fluctuation 
data appears to be that MB dijets play a dominant role 
in all high-energy nuclear collisions, consistent with pre¬ 
vious analysis of yields, spectra and correlations at the 
RHIC. Jet manifestations are clearly evident in extensive 
measures n c h and Pt and are simply and accurately rep¬ 
resented by the TCM over a range of collision systems 
and energies. Measures that rely on ratios of means such 
as pt or ratios of random variables such as (p t ) present 
an ambiguous picture because of possible cancellation of 
dijet contributions. Data in the form \/~C/pt, a ratio of 
ratios including an additional factor 1 /n 2 ch compared to 
conventional variance fluctuation measures, are difficult 
to interpret in that form. Transformation of such data to 
a variance-difference format presents a clearer picture. 

This material is based upon work supported by the 
U.S. Department of Energy Office of Science, Office 
of Nuclear Physics under Award Number DE-FG02- 
97ER41020. 


Appendix A: pt TCM for A-A collisions 

The TCM for A-A collisions is based on the Glauber 
model in which the fractional cross section (centrality) 
ct / cq is related to geometry parameters N part the num¬ 
ber of projectile nucleon participants, Nf nn the number of 
binary N-N encounters and v = 2N part /Ni, in the mean 
participant pathlength in number of N-N encounters. 
The correspondence with observable n c h can be estab¬ 
lished from the MB cross-section distribution on n c h- 
For the present study the correspondence between n c h 
reported in Ref. [T] and Glauber model parameters was 
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determined as described in Sec. IV C| and Ref 


For A-A collisions the pt TCM of Eq. (|23|) or (24) 


must be modified in three ways: (a) the multiplicity hard 
component increases with centrality as nu(v), (b) due to 
modified parton fragmentation in more-central A-A colli¬ 
sions the spectrum hard-component shape changes (soft¬ 
ens) with centrality leading to variation of p t} h as Pt,h{ v ) 
and (c) an N-N “first encounter’ effect must be accom¬ 


modated, with details presented in Sec. IV B 


The direct extensi on o f p-p n c h = n s + rih to A-A is 
the first line of Eq. (All where the N-N soft and hard 


components are sealed up by the corresponding Glauber 
parameters. The observed trend for hadron production 
in A-A collisions implies that nh for the first N-N en¬ 
counter is the same as that for p-p independent of the 
A-A centrality, but for v — 1 subsequent encounters nh 
transitions to n' h (y) depending on A-A centrality. The 
consequence is the second line that accurately describes 
n c h trends for a variety of collision systems 


N 


part 


Tick — HsiNpart/ty ~\~ ^hiy^-^bv 

-n ch = n pp [l + x(v)(v — 1)], 


(Al) 


hh{i / ) is an average over all v N-N encounters whereas 
n' h (y) or x{y) applies only to the v — 1 subsequent or 
secondary encounters. 

For a self-consistent description the same argument 
should be applied to Pt,h{y) such that in the first N-N 
encounter the p-p value holds while thereafter the value 
may change. The A-A TCM for pt{v) with p t cut is then 


Pt,AA( v ) 


^sPt.,3import/‘^‘) T" 'fi'h{y') , Pt,h{y)Nbin 

N s (Np ar t / 2 ) -)- flh{y^}Nhin 
Pt,pp + x{v)p t , h {v){y - f) 

ri pp /n pp +x(v) (v - 1 ) 


(A2) 


where n' pp /n pp ~ 0.75 for p tiCU t ~ 0.18 GeV/c. For 
a complete A-A p t TCM description it remains to de¬ 
fine quantities x(v) and p ti h{v). x{y) is defined in the 
second line of Eq. (22) with uq = 2, xq = 0.028 and 
x\ = 0.141 for 2.76 TeV Pb-Pb collisions. The hard 
component pt,h(v) for the same system is defined by 


p t h {v) = im + 1.70{1 - tanh[(zz - i^i)/0.42]}/2(A3) 


where n pp = n s + nh and x(v) = n' h (v)/n pp . Note that with zq = 1.75, all as reported in Ref. [13] . 
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